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Abstract. Results concerning the statists of r-adic processes and their fractal 
properties are reviewed. The connection between singular eigenstates of the 
statistical evolution of such processes and popular fractal curves is emphasized. 

In memory of our friend and colleague Shuichi Tasaki whose untimely death cut short 
a remarkable scientific career 

Among the many important scientific achievements of Shuichi Tasaki were 
his contributions to the statistical properties of piecewise linear maps and the 
characterization of their eigenstates. In particular his work was key to understanding 
the role played by singular measures in the statistical evolution of chaotic maps and 
its connection to the mathematics of fractals [U [2]. This field, which attracted 
much attention in the non-equilibrium statistical physics community in the mid- 
1990's [2111], had been popular among many Japanese mathematicians, in particular 
early in Tasaki's scientific career [5]. In this respect, it is perhaps not surprising 
that one of Tasaki's favorite examples of fractals, which helped understand how the 
non-equilibrium states of volume-preserving strongly chaotic systems acquire fractal 
properties [6], was actually introduced by Teiji Takagi 7j, the founder of the school 
of modern mathematics in Japan, who, more than a century ago, had proposed it as 
a simple example of a continuous but nowhere differentiable function. 

It is the purpose of this article to review some of Tasaki's contributions to this 
field, and draw observations which aim to underline the similarities between the so- 
called hydrodynamic modes of diffusion of simple model systems and some popular 
fractals, which include the von Koch and Levy curves. 

1. Singular measures of dyadic processes 

The statistical properties of coin tosses are full of mathematical wonders and Shuichi 
Tasaki was well aware of it. Thus assume a fair coin tossing game and a sequence 
{uji, . . . ,u k } of k binary digits, £ {0, 1}, i = 1, . . . , k, where say the symbol 
stands for heads and 1 for tails. The probability measure which assigns probability 
1/2 to every symbol irrespective of which of heads or tails came at the previous toss is 
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invariant under such a process, which in particular is to say that every set of k binary 
digits has the same measure 1/2 . No surprise there. 

Consider however an arbitrary number p, < p < 1, and let /j, p denote the 
probability measure which assigns measure p to heads and 1 — p to tails. This measure 
is itself invariant under fair coin tossings [S] . 

The reason is simply that the set . . . , w&} is the union of the two disjoint sets 
of length k + 1, 



Since /x p ({0, oji, . . . , oj k }) = pn P {{u>i, . 



,w fc }. 

. ,u)k}) and /j, p ({l,U)i,. 



(1.1) 

■,w fc }) = (1 - 



p)fj,p({uJi, . . . , u>k}), we have that the probability measure of the left and right hand 
sides of equation ( 1.1 ) is the same, irrespective of the choice of the value p. We might 



say, in other words, that fj, p is as good an invariant measure as the uniform measure 
Mi/2 is- 

There is however an essential difference here, which is that H1/2 is the natural 
invariant measure for the fair coin tossings. For p ^ 1/2, fi p is in fact a singular 
measure, which can be thought of as a singular Lebesgue function [3]. Using the 
isomorphism between fair coin tossings and the angle-doubling map, x H ► 2x(mod 1), 
it is a simple exercise to identify fi p with the measure lifted on the corresponding 
cylinder sets of the unit interval, whose cumulant fi p ([0,x]) is the unique function 
f p (x) satisfying 

pf p (2x), 0<x< 1/2, 

(l-p)f p (2x-l)+p, 1/2<x<1. 

Except for p — 1/2, for which fi/2(x) = x, f p is strictly increasing and continuous, 
but has zero derivatives almost everywhere with respect to the Lebesgue measure |10) . 



f P (x) 



(1.2) 



A specific example is shown in figure 1(a) 
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Figure 1. (a) Lebesgue's singular function f p , equation (1.2 I, and (b) fp 



(l-p) 



equation | |1.3| , both with p = 
curves are computed over 2 14 



2/3. Here and below, 
points. 



unless otherwise stated, the 



Notice however that [i p is in fact the natural invariant measure of the biased coin 
tossing which gives probability p to heads and 1 — p to tails. Its cumulant corresponds 
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to a uniform measure, which can be identified as the solution of the more general 
functional equation 



fi q \x) 



P/^(f), 

(l-p)/i 9) (f5f)+P, 



< x < q, 
q < x < 1. 



(1.3) 



Setting q = p indeed yields fp P \x) 



On the other hand, the function fp 



(1-p) 



illustrated in figure 1(b) is of the Lebesgue singular type and arises as the cumulant 
of the natural invariant measure of a dissipative baker map, projected along the 
contracting direction [9]. 



2. Complex- valued measures of dyadic processes 



Though equation (1.2 1 is a functional equation characterizing a probability measure, 



it is not in itself restricted to real values of the parameter p. The uniqueness of its 
solutions is indeed warranted for every complex parameter p such that \p\, |1 — p\ < 1 
[TTj [12]. In such cases the solutions of equation (1.2 1 are self-similar sets on the 



complex plane and it is a simple calculation to show that their Hausdorff dimension 
is the solution c?h of the following equation [13] , 

(2.1) 



\ P \ d * + \i - P \ a » = i. 

A graphical representation of its solution is shown in figure [2] 



,l<fa _ 



d„(p) 




Figure 2. Hausdorff dimension l |2.1[ l of the curves f p defined through equation 

JTJ}. 

Note in particular that the parameter values p — 1/2(1 + e lv ), < ip < 2ir, 



define a circle of complex parameters such that the dimension (2.1 ) is equal to djj = 2 



Curves with parameters whose values lie inside that circle have Hausdorff dimension 
between 1 and 2. Curves with parameters which verify the conditions \p\ < 1 and 
|1 — p\ < 1, but which lie outside that circle have Hausdorff dimension 2 [12] 
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As shown by Tasaki and collaborators [2], the functions f p with complex 
parameter 

p{k) = = i + t tank), (2.2) 

2 cos fc I 

with fc a real number such that |fc| < vr/3, occur as representations of the left 
eigenvectors of the Frobenius-Perron operator of an expanding piecewise linear map 
of the real line related to the angle-doubling map. 

The same curves also occur as special cases of hydrodynamic modes of diffusion, 
associated with multi-baker maps on a ring |14j . where k is the associated wavenumber. 
The dynamics is defined according to 

/ \ J (n-l,2z,§), 0<x< 1/2, 

(„ + l,2z-l,^i), 1/2<z<1, (2 ' 3) 

where n = 0, . . . , N — 1 and n ± 1 are understood to be taken modulo AT. Given an 
initial distribution, the relaxation of statistical ensembles to the uniform equilibrium 
measure takes place exponentially fast at rate ~ 7V~ 2 and is best characterized in terms 
of the cumulant measure ^t(n, [0, 1], [0,y]) of phase points (n, x, y), with < x < 1 
and < y < y, after t iterations: 

m(n, [0, 1], [0, y}) = cos t (fe)F fc (y)e tfc ", (2.4) 

where fc = 2i:m/N ', m = 0, . . . , AT — 1 are the wavenumbers, the coefficients are set 



by the initial distribution, and are solutions of the system of equations (1.2) with 



the choice of parameter (2.2), Fk(x) = f p (k){x), viz. 



,F k (2x), 0<x<l/2, 



I fef^(2«-l) + 5£l. 1/2<«<1. 
The dimension ( |2.1[ ) of these curves can be computed explicitly [2] : 

d H (fe) = ^ g2 (2.6) 
m ; log(2|cosfc|) v ; 

In particular, for fc = 7r/4, we obtain a curve of Hausdorff dimension dn(n /4) = 2 
which covers positive areas of the plane and corresponds to the Levy dragon [151 . 
Another remarkable value of the parameter is fc = it/ 6, for which we have the 
dimension d H ( 7r /6) = log 4/ log 3. These two curves are displayed in figure [3] 

In the context of hydrodynamic modes of diffusion |16j , the complex curves defined 



by equation (1.2| acquire a clear physical meaning , which comes about when the 



parameter fc in (2.2) is small: |fc| 1. The Hausdorff dimension (2.6) then becomes 



which provides a relation between the Hausdorff dimension of the hydrodynamic 
modes, the diffusion coefficient of the process and the positive Lyapunov exponent 



of the underlying dynamics (2.3). Furthermore, considering the functional equation 



(2.5), we obtain the equivalent of a gradient expansion for F k |17j . 

F k (x) = x + ikT(x) +C(fc 2 ), (2. 
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(b) Levy dragon with dimension 2 



Figure 3. Specific examples of hydrodynamic modes for the dyadic process 
defined by the functional equations j2.5| ■ 



where T{x) is the Takagi function [7], which can be defined through the functional 
equation: 

< x < 1/2, 



See figure [4] 



T(x) 



x + \T{2x) 



1 - x + \T(2x - 1), 1/2<x<1. 



(2.9) 
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Figure 4. The graph of the function F^, equation | |2.5| , with k = 1/10 (a), 
approaches that of the Takagi function ||2.9| (b). 



The result (2.8 1 should be compared with that of Hata and Yamaguti [TBI I19j. 
who recover the Takagi function as the derivative of the curves f p , equation ( 1.2 ), with 
respect to their parameter, T(x) = df p {x)/dp\ p =i/2- 



3. Levy, von Koch and Heighway curves 

The Levy dragon shown in figure [3(b)] is but one example in a long list of popular fractal 
curves which are most often defined recursively through iterated function systems 
[2"U] or using L-systems [5T]. The von Koch curve [22| and the Heighway dragon, 
popularized in Martin Gardner's Mathematical games [23], are two other similar 
examples. In analogy with equation (1.2 1, it is simple to identify linear contractions 
which generate these sets [2H [23 ■ 
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The solutions of equation (1.2) with complex parameters produce a family of 
curves somewhat similar to the Levy dragon. Looking at the iterative construction of 
these solutions, we start from the initial line segment joining Fk(0) = and -Pfc(l) = 1 
and obtain i*fe(l/2) = e lfc /(2cosfc). At the next iteration, we add the two points 
Ffc(l/4) and Ffc(3/4). The result on the complex plane is a collection of two isosceles 
triangles whose long edges correspond to the two smaller edges of the triangle formed 
by the three initial points -Ffe(O), Fk(l/2) and -Ffc(l). The next iteration produces four 
smaller isosceles triangles which are stacked upon the small edges of the two existing 
ones. See figure 5(a)| 




(a) 



(b) 



Figure 5. Comparison between the first few steps of the iterative constructions 
of (a) the Levy curve and (b) the von Koch curve. 

A similar construction, which consists of stacking up the triangles alternatively 
upwards and downwards from one iteration to the next, as shown in figure |5(b)| 
produces curves similar to the von Koch curve. They are most easily obtained as the 



self-similar sets associated with a linear contraction similar to equation ( 1.2 ), but that 
further involves complex conjugation: 



P, 



0<x< 1/2, 
1/2 < x < 1, 



(3.1) 



(1 - p)g* p (2x - 1) 

where * denotes the complex conjugation. Equations (1.2) and (3.1) share the same 



real solutions but have different sets of solutions for complex parameter values. 



Nonetheless they share the same Hausdorff dimension (2.1). As pointed out by de 
Rham [T2], the choice p = 1/2(1 + i/yS) produces the von Koch curve. 

Taking the parameter p as defined by equation (2.2), we write Gk{x) = g p (k)( x )> 
for which we have the functional equation: 



G k (x) 



2 cos k 



2 cos k 



Gt(2x), 
G* k (2x-1) + 



(3.2) 



2 cos k) 



< x < 1/2, 
1/2 < x < 1. 

7r/6 corresponds to the von 



As illustrated in figure |6j the choice of parameter k -- 
Koch curve. For k = tt/4 we obtain a Peano-like curve. 

As such, the functional equation (3.2) is different from equation (2.5) so its 



solutions do not correspond to hydrodynamic modes of the dyadic multi-baker map 
(2.3). It is however straightforward to see that they are the hydrodynamic modes 
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Figure 6. Specific examples of solutions of the functional equation | |3.2| l yield 
(a) the von Koch curve and (b) a Peano-type curve. 



associated with the four-adic map, 



(n,4x, |), 



(n,x,y) i-» 



1,4a: -1,^), 
1,4a: -2,^ 

3, 



< x < 1/4, 
1/4 < x < 1/2, 
1/2 < x < 3/4, 
3/4 < x < 1, 



(3.3) 



(n- 
(n- 

(n,4ar- ~, 4 

with wavenumbers halved. Note that the iterative construction of the von Koch curve 
using a four-adic process amounts to skipping every odd step in the dyadic-based 
iterative construction shown in figure 5(b) The former is often preferred over the 



latter in the literature, see for instance [26]. even though the dyadic representation 
based on equation (3.2) is indeed the most compact. 



The Heighway dragon and related curves can also be obtained through functional 



equations similar to equation ( 1.2 ). They are based upon iterative constructions which 



combine both upwards and downwards triangles, as shown in figure [7] 




(a) 





(b) 



Figure 7. First few steps of the iterative constructions of asymmetric curves 
similar to the the Levy and von Koch constructions of figure [5] 

The Heighway dragon is in fact a particular case of a hydrodynamic mode of a 



dyadic map similar to the multibaker map (2.3), where the angle-doubling map is 
replaced by the tent map x h- > 2x if < x < 1/2 and 2(1 — x) if 1/2 < x < 1, namely 

(71-1,1-201,1), 0<x<l/2, 
(n+ 1,2a: -1,1- |), 1/2<x<1. 



(n,x,y) i-> 



(3.4) 
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Note that the tent map here appears upside down along the expanding direction. 
This choice of combination of maps along the x and y coordinates is taken so 
the map has the property of being time reversal symmetric under the induction 

(n,x,y) i ^ (n, 1 



which is also the case of the multi-baker map (2.3) 



Another such map, very similar to equation (|3.4|), is 
(n,x,y) i y 



(3.5) 



(n-l,2x,±=^), 0<x<l/2, 
(n + l,2(l-a;),i!±i), 1/2 < x < 1. 

Here the usual tent map acts along the expanding direction. For the sake of our 



argument however, we prefer using (3.4) since the tent map appears along the 



contracting direction under the time-evolution of phase-space densities. 



Using, for the evolution of statistical ensembles under (3.4), an expansion similar 



to equation (2.4), we identify the modes 



P, 



2 cos k 



- : Pk(2x), 



1 - 



2 cos k 



P k (2(l-x)), 



< x < 1/2, 
1/2 < x < 1. 



(3.6) 



The value k = ir/4 yields the Heighway dragon, as shown in figure |J The map (13. 5| 
produces modes which can be obtained from equation (3.6) by a simple symmetry. 
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Figure 8. Examples of different solutions of equation ( |3.6[ l with (a) k = tt/6 and 
(b) k = 7r/4. The latter is knwon as the Heighway dragon 



Finally, by analogy with equation (3.2), we obtain another set of associated 



fractals by complex conjugating the function on the right-hand side of the functional 
equation, 



2 cos fe 



Q%(2x), 



< x < 1/2, 



2 cos k 



£(2(1 -a)), 1/2<x<1. 



(3.7) 



Examples of solutions are shown in figure [9j Here again these curves can be identified 
as the hydrodynamic modes of a four-adic map. 



4. Complex measures of r-adic processes 



The above considerations are not limited to dyadic processes. Consider for instance 
the iterated systems shown in figure [10] The functional equations whose solutions 
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Figure 10. These three examples of iterative constructions of fractal curves are 
based upon triadic processes which can be conveniently written under the form 
of functional equation systems ||4,l[l-||4~3|l. 
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reproduce these familiar curves are based upon triadic processes: 



17* (3s) 



U k (x) 



V k (x) 



W k {x) = 



1+2 cos k 



;U k (3x-l), 
: U k (3x-2), 



1 + 2 COS k 

e 

1+2 cos k 
e lfc + l 

1+2 cos k 1 1+2 cos k ' 
l+2cosfc^4*( 3x )' 

iq$^*(3*-2), 



1+2 cos k 

e lfc + l 
1+2 cos k 



1+2 cos k 

e' fc 
1+2 cos fc 
e' fc + l 



W k (3x), 



1+2^*^(3* 
I^kW k (3x- 



'1) 
2). 



< x < 1/3, 

1/3 < x < 2/3, (4.1) 

2/3 < x < 1, 

< x < 1/3, 

1/3 < x < 2/3, (4.2) 

2/3 < x < 1, 

< x < 1/3, 
1/3 < x < 2/3, (4.3) 
2/3 < x < 1, 



1+2 cos k 

defines the hydrodynamic modes of diffusion of the triadic multi- 
. Equations (14. 21) and (4.3 1 correspond, on the other hand, to the 



Equation ( |4.1[ 
baker map [14 

hydrodynamic modes of nine-adic multi-baker maps associated with random walks 
with assigned probabilities 1/9 to jump by two units to the left or right, 2/9 by a 



single unit, and 1/3 to remain put. Notice in particular that equation (4.2) produces 
the Sierpinski gasket 



for k — 7r/3, as shown in figure 11 




Fig ure 11. Examples of different solutions of equations (a) | |4.1[ ), (b) j4.2| , and 
(c) fljjb, with k = 7r/3, here computed over 3 9 points. 



Further examples of piecewise linear maps and their associated eigenstates can 
be found in 1251 . 



5. Concluding remarks 

In [6], Tasaki and Gaspard showed that the Takagi function describes the non- 
equilibrium stationary state of a multi-baker map, and used its properties to show 
that this system obeys Fick's law, which provides an expression of the mass current 
as minus the product of the mass density gradient and diffusion coefficient. This 
result thus provided a derivation of a phenomenological law of thermodynamics in 
terms of the phase-space dynamics of the model, which subsequently lead Gaspard to 
identify the fractality of the non-equilibrium stationary states as the source of entropy 
production [21] . A recent survey of these results with further extensions can be found 
in [SD]. 

In this paper we have presented some of these results in the more general 
context of fractal curves and self-similar sets defined by linear contractions, and 
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emphasized their connections with the statistical properties of maps that generate 
the associated functional equations. Shuichi Tasaki contributed in important ways 
to our understanding the connection between singular functions and the statistical 
properties of r-adic maps. We are deeply saddened that we are not able to talk to 
Shuichi about the excursion into the world of fractals described here. We know that 
we would have learned still more about them, and benefited from his careful and clear 
explanations. 
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